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Three-body nuclear reactions in two-nucleon plus core systems are described in the framework 
of exact scattering equations including the core excitation. A nucleon-core optical potential is 
constructed that can be easily adjusted to the reference potential and thereby to the experimental 
two-body data, if available. This constitutes an important improvement over the simple deformation 
of the potential used previously that violated the original fit to the data. Predictions for elastic, 
inelastic, and transfer reactions involving ^°Be and ^“^Mg nuclear cores are obtained. The new optical 
potential leads to a moderate increase of cross sections. 

PACS numbers: 24.10.-i, 21.45.-v, 25.45.Hi, 25.40.Hs 


I. INTRODUCTION 

Deuteron {d) scattering from a nucleus A, consisting 
of A nucleons, microscopically is an A -f 2 body prob¬ 
lem. While various A > 3 cases have been considered 
using approximate methods W , a rigorous solution of ex¬ 
act scattering equations p-Q has been achieved so far 
for three- and four-body systems only [m. In addi¬ 
tion, there exist also numerous approximate treatments 
for three-body reactions [iMl- To apply the avail¬ 
able three-body techniques to deuteron-nucleus scatter¬ 
ing, this process is often approximated by a three-body 
problem where the nuclear core A is treated as a struc¬ 
tureless particle whose interaction with the proton (p) 
and the neutron (n) is given by complex or real poten¬ 
tials. In some cases this may be a reasonable approach, 
but in others it is necessary to go beyond the picture of an 
inert core A and take into account its internal degrees of 
freedom. It has been shown in several works |16l - [l^ that 
the excitation of the core may be an important reaction 
mechanism and needs to be taken into account. There 
are several formulations of rigorous three-body scattering 
equations including core excitation (CeX) [1^ , how¬ 

ever, numerical results were obtained only in Ref. M- 
As dynamic input these calculations use nucleon-nucleon 
(NN) and nucleon-core (A^A) interactions. There is a 
number of standard parametrizations for NA optical po¬ 
tentials (OP) without CeX. The CeX is usually included 
by deforming these standard potentials to allow the cou¬ 
pling between ground (A) and excited (A*) states of the 
core. However, this way the additional TV A* component 
in the two-body scattering equation together with the 
potential deformation distort the elastic NA amplitude 
that deviates from the original elastic amplitude calcu¬ 
lated using standard potential without CeX. The latter 
is usually fitted to the A{N, N)A elastic experimental 
data, thus, the description of the data gets destroyed 
when including the CeX via the simple deformation of the 
potential. Furthermore, if the NA — NA* coupling po¬ 
tential was fitted to the inelastic A{N, N')A* data using 
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the distorted-wave Born approximation (DWBA) which 
is the usual case, the fit to A{N,N')A* data is lost as 
well. Thus, an additional adjustment of the OP param¬ 
eters needs to be performed to restore the description of 
TV A elastic and inelastic scattering data. This implies a 
new fit procedure for each nucleus at each energy, which 
may be quite tedious. On the other hand, one may take 
the advantage of already existing fits to the experimen¬ 
tal data obtained using DWBA and standard OP. For 
this a TVA potential with CeX is needed that, when in¬ 
serted into the exact two-body multichannel Lippmann- 
Schwinger equation, yields the standard OP results for 
the A(TV, TV)A elastic scattering and DWBA results for 
the A(TV, TV')A* reaction. In the present work I propose 
a method for achieving this goal, i.e., I construct a TVA 
potential with CeX that at a given energy exactly re¬ 
produces the standard OP and DWBA amplitudes for 
A(TV, TV)A and A(TV, TV')A* scattering processes, respec¬ 
tively. Thus, the desired consistency with the TVA ex¬ 
perimental data is ensured. In case there are no data 
available, this approach still allows for a more precise 
evaluation of the CeX effect, not affected by the OP mis¬ 
match. In three-body systems I aim to study the changes 
in observables caused by this improvement of the OP. I 
therefore reanalyze the reactions calculated in Ref. [l^ 
using the simple deformation of the OP as well as present 
several new cases. 

In Sec. mi derive the TVA potential with CeX, and in 
Sec.imil recall the three-body scattering theory. Exam¬ 
ple results for d + ^°Be, p + ^^Be, and d + '^^Mg reactions 
are presented in Sec. IIYI while the summary is given in 
Sec. El 


II. SUBTRACTION METHOD FOR 
NUCLEON-CORE POTENTIAL 

I employ the extended Hilbert space TLg © Hx where 
the two sectors correspond to the core being in its ground 
(g) or excited (x) state [l^. I aim to construct the TVA 
potential coupling the two sectors and denote its compo¬ 
nents by Vgg, Vxx, and Vxg = they are represented 
graphically in Fig. [TJ The coupled-channel Lippmann- 
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FIG. 1. Nucleon-core potential in the three-body Hilbert 
space. The particles are represented by vertical lines, the 
core in its excited state is distinguished by thick lines, the 
potential by horizontal dashed lines. The diagrams (a), (b), 
(c), and (d) correspond to Vgg, Kg, Vgx, and V^x- 

Schwinger equation 

('^g 9 '^gx\ _ f^gg^gx\ , f^gg^gx\ q f'^gg'^gx\ 
[Txg Txj [Vxg Vxj ^ Vxj ° [Txg T, ^ 

( 1 ) 

yields the respective components Tji of the NA tran¬ 
sition matrix at the available energy E with Gg = 
{E + iO — Hq)~^ being the free resolvent. I emphasize 
that the extended free Hamiltonian Hq besides the kinetic 
energy operator contains also the internal core Hamilto¬ 
nian whose contribution vanishes for Tig but is equal to 
the core excitation energy (toa* — ttia) for T-Lx- 

At energy E = Eg the component Tgg describing the 
NA elastic scattering is demanded to reproduce the tran¬ 
sition matrix 


FIG. 2. Diagrammatic representation of the contributions 
at the two lowest orders subtracted from the nucleon-core 
potential. 

where Vi are the respective projectors and 82 = l 32 Ro is 
the deformation length. 

To fulfill the above demands, Txg in Eq. o is resolved 
as 

Txg = (1 - VxxGo)-^Vxgil + GoTgg) (5) 

and used in this for the Tgg component, leading to 

Tgg = Vgg{l + G0Tgg) + VgxG0{l-VxxG0)-Wxg{l + G0Tgg). 

( 6 ) 

Furthermore, I make use of the identity 

— VxGq) ^ = 1 + txGo (7) 

obtained from Eq. ©. By comparing Eqs. ([6]) and ([5]) 
with @ and (|3]) one observes that the conditions Tgg = tg 
and Txg = T^^^^ are satisfied at E = Eg by choosing 


tg=VgE VgGotg (2) 

obtained with the NA potential Vg not including CeX, 
usually taken from one of the standard parametrizations. 
There is no such demand for Txx due to the lack of the 
experimental data for the NA* elastic scattering. The 
respective transition matrix for the core in its excited 
state but without coupling to the ground state 

tx — '^x T '^xG^tx (3) 

is obtained with the NA* potential Vx that is not con¬ 
strained by the data but is usually taken from the same 
standard OP parametrization as Vg. 

Furthermore, at E = Eg I also require Txg to reproduce 
the DWBA amplitude 

= (1 + GGo)E,™(l + Gotg). (4) 

The potential coupling the two Hilbert sectors 

Jig and T-Lx usually is obtained by deforming the cen¬ 
tral part of Vg. In the rotational model [2l| one as¬ 
sumes that the core has a permanent quadrupole de¬ 
formation and replaces in Vg the nuclear radius Rq by 
R = i?o(l + /^ 2 d 2 o(C)); where (32 is the quadrupole de¬ 
formation parameter and ^ describes the internal core 
degrees of freedom in the body-fixed frame. If the cen¬ 
tral part of Vg is a function of (r — Rq), e.g., the Woods- 
Saxon function, is given by VxVg{r — 62 Y 2 o{i))Vg 


Vgi 

Vx, 


Vx, 


~ ^g VgxQ^i^ txgQ)Vxg, 

_ wDWBA 
^xg J 

= Vx, 


(8a) 

(8b) 

(8c) 


with 5 q = {Eg -I- fO — Hq)~^. Thus, the essential idea of 
the present method is subtracting from the elastic ampli¬ 
tude the contributions that are explicitly generated in the 
scattering equations by the coupling to the core excited 
state. For this reason it is called the subtraction method. 
The subtracted contributions at the two lowest orders are 
diagrammatically represented in Fig. [2] These contribu¬ 
tions are nonlocal in the coordinate space but this has 
no disadvantage when the calculations are performed in 
the momentum space. 

Note that a similar subtraction method has been used 
in the past to readjust the purely nucleonic part of the 
two-nucleon potential with explicit A-isobar excitation 
[ 2 ^ , thereby improving the fit of the elastic two-nucleon 
scattering data. The present method is more general in 
the sense that it also hxes the inelastic amplitude. 


III. THREE-PARTICLE SCATTERING 
EQUATIONS 

Exact three-body scattering equations of Faddeev or 
Alt-Grassberger-Sandhas (AGS) type including the GeX 
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have been discussed in Refs. [a I 20 II : their practical so¬ 
lution was implemented in Ref. [l9|. In the extended 
Hilbert space they acquire the standard form of the AGS 
integral equations for three-body transition operators 

3 

Ufia = ^jia Gg ^ -f Tj G^U-fa (9) 

7=1 

that couple the two sectors T-Lg and T-Lx much like the 
two-body transition operators in Eq. ([T]). The sub¬ 
scripts a, /3 ,7 label the spectator particles (interacting 
pairs in the odd-man-out notation) for operator com¬ 
ponents, while = 1 — 5Pa- Following the develop¬ 
ments of Ref. [a|) I solve Eq. m numerically in the 
momentum-space partial-wave framework, including the 
proton-core Coulomb force via the screening and renor¬ 
malization method [23l - l^ . The scattering amplitudes 
are given by the on-shell matrix elements of the transi¬ 
tion operators Upa calculated between initial and final 
channel states. 


IV. RESULTS 

In this section I present results for two-cluster elas¬ 
tic, inelastic, and transfer reactions initiated by p-|-^°Be, 
d ^°Be, p ^^Be, and d -I- ^“^Mg collisions. For both 
^°Be and ^^Mg cores the ground (first excited) states 
have spin and parity O'*" (2+), while the respective ex¬ 
citation energies are 3.368 and 1.369 MeV. The main 
goal is to evaluate the changes in the observables due 
to the use of the improved nucleon-core potential includ¬ 
ing the CeX with subtraction. I present three types of 
calculations. All of them use a realistic CD Bonn po¬ 
tential [ 2 ^ for the np pair, but differ in NA interac¬ 
tions: (i) single-particle (SP) model for the core, i.e., 
neglecting CeX in NA potentials, labeled as SP in the 
following; (ii) including CeX via simple deformation of 
the SP potential without subtraction and readjustment, 
labeled CX(no subtr) in the following; and (iii) includ¬ 
ing CeX using improved NA potentials with subtraction 
and proper readjustment as described in Sec. |lTl labeled 
CX in the following. Thus, the difference between CX 
and SP will yield the the CeX effect, while the difference 
between CX and CX (no subtr) will evaluate the im¬ 
portance of the potential improvement. The calculations 
SP and CX(no subtr), both based on the Chapel Hill 89 
(CH89) parametrization (2^, are the same as in Ref, [j^. 
but CX(no subtr) was labeled simply CX in Ref. [T^. I 
therefore only have to describe the calculations CX with 
subtraction. For the nA pair in the partial waves with 
{An) bound states I take over the real binding poten¬ 
tials from Ref. [l^ since they are already adjusted to 
experimental binding energies and need no subtraction. 
Depending on the partial wave, these potentials may sup¬ 
port deeply bound states | 6 o) that are Pauli forbidden 
and therefore have to be projected out; this is achieved 
by adding a separable term |&o)P(|&o| with P > 1 GeV to 



FIG. 3. (Color online) Differential cross section for 
^°Be(p,p)^°Be and ^°Be(p,p')^°Be* reactions at Ep = 39 
MeV. Results including the CeX with and without subtrac¬ 
tion, i.e., CX and CX(no subtr), are given by solid and 
dashed-dotted curves, respectively. CX predictions coincide 
with SP results neglecting the CeX for the elastic scattering 
and with DWBA results for the inelastic reaction. The elastic 
experimental data are from Ref. [^ . 


the local nA potential and thereby movin g th e Pauli for¬ 
bidden state 1 60 ) to a large positive energy [^. For nA in 
other partial waves and for pA in all partial waves I take 
the potentials including CeX with subtraction; they are 
derived from the CH89 parametrization unless explicitly 
stated otherwise. The subtraction energy Eg coincides 
with the energy at which the OP is taken, i.e., half of 
the deuteron energy for nA but proton energy Ep for 
pA. However, if the nA potential is complex in all partial 
waves as in the case of the present d-|-^"'^Mg calculations, 
the p + (An) channel is absent and therefore the pA po¬ 
tential is taken at Ed/2 as well. These choices are the 
same as in Ref. M- The employed quadrupole deforma¬ 
tion parameters for V-^°Be interactions are ,82 = 0.67 
and 62 = 1.664 fm as in Ref. [l^, while for V-^^Mg I 
use (32 = 0.5 and 62 = 1.69 fm. In a few selected cases I 
present CX results derived from Watson [ 2 ^ and Koning- 
Delaroche parametrizations; those calculations use 
the S2 values listed above and are labeled as CX(W) and 
CX(KD), respectively. 

First in Fig. [3] I demonstrate the importance of the 
proper potential readjustment by comparing the CX and 
CX(no subtr) results in the two-body system. As an 
example I show the differential cross section da/dfl for 
elastic and inelastic p + ^°Be scattering at Ep = 39 MeV 
proton energy as a function of the center-of-mass (c.m.) 
scattering angle 0c.m.- Here and in the following the 
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FIG. 4. (Color online) Differential cross section for p+ ^^Be 
elastic scattering at Ep = 35.3 MeV. Results of SP, CX, and 
CX(no subtr) potential models are given by dotted, solid, and 
dashed-dotted curves, respectively. The experimental data 
are from Ref. [^ . 



FIG. 5. (Color online) Differential cross section for transfer 
reactions ^°Be(d,p)^'^Be at Ed = 21.4 MeV leading to the 
ground (|^) and excited (^ ) states of ^^Be. Curves are as 
in Fig. |4] and experimental data are from Ref. [s^. 


elastic differential cross section is given as a ratio to the 
Rutherford cross section dan/dVl. The differences be¬ 
tween the CX and CX(no subtr) calculations are of mod¬ 
erate size for elastic scattering but become more signif¬ 
icant for the ^°Be(p,p')^°Be* reaction, reaching nearly 
40% at forward angles. I remind that by construction 
the CX and SP results coincide for the elastic scattering 
while the CX and DWBA results coincide for the inelastic 
scattering. 

Observables for p + ^°Be and p + ^^Be elastic scatter¬ 
ing are quite strongly correlated. Consistently with this 
fact the SP and CX models, yielding identical results for 
^°Be(p,p)^°Be, agree quite well also for ^^Be(p,p)^^Be 
at Ep = 35.3 MeV as shown in Fig. 01 indicating that 
the CeX effect on the p -I- ^^Be elastic cross section is 
very small. However, the p-^°Be potential, that does not 
properly describes elastic two-body data as in the case of 
CX(no subtr), overestimates the CeX effect. 
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FIG. 6. (Color online) Differential cross section for 
^^Be(p,d)^°Be transfer reactions at Ep = 35.3 MeV leading 
to the ground (O'*") and excited (g”*") states of '^'^Be. Curves 
are as in Fig. [4] and experimental data are from Ref. [^ . 


Next I show in Fig. [5] the differential cross sections for 
the ^°Be(d,p)^^Be transfer reactions at Ed = 21.4 MeV 
deuteron energy. Quite surprisingly, the improvement 
of the 7V-^°Be potentials has a very small effect on the 
transfer to the ^^Be ground state i"*". For the transfer 
to the ^^Be excited state ^ the effect is of moderate 
size but not beneficial. Thus, the ^°Be((i,p)^^Be* reac¬ 
tion still remains an unresolved problem and calls for a 
more sophisticated model. In Fig. [6] I present results for 
^^Be(p, c?)^°Be transfer reactions at Ep = 35.3 MeV; it 
corresponds to ^°Be((i,p)^^Be at Ed = 40.3 MeV. At this 
higher energy the importance of the proper potential ad¬ 
justment is more visible, especially for the reaction lead¬ 
ing to the 2+ excited state of the ^°Be core. The latter 
finding is not unexpected given the results for p -h ^°Be 
inelastic scattering in Fig. [31 For both ground and ex¬ 
cited ^°Be states the description of the experimental data 
is improved, although some discrepancies still remain. In 
fact, these predictions are quite sensitive to the choice 
of the optical potential from which the CX model is de¬ 
rived. I illustrate this finding in Fig. |3 by comparing 
the CX results for the ^^Be(p, (i)^°Be differential cross 
sections based on CH89, Watson, and Koning-Delaroche 
OP parametrizations. I admit that no one of them re¬ 
produces the experimental data for both reactions simul¬ 
taneously: transfer to the ^°Be ground state 0+ is best 
described by CX(KD) while to the excited state 2+ by 
CX(W). 

In Ref. [l^ it was demonstrated that the CeX effect in 
transfer reactions is, in general, much more complicated 
than just a simple rescaling of the SP differential cross 
section by the respective spectroscopic factor (SF). This 















5 



0 30 60 90 

®c.m. (deg) 


FIG. 7. (Color online) Differential cross section for 
^^Be(p,d)^°Be transfer reactions at Ep = 35.3 MeV lead¬ 
ing to the ground (0”*") and excited (2^) states of ^'^Be. CX 
predictions based on CH89 (solid curves), Koning-Delaroche 
(double-dotted-dashed curves), and Watson (dashed curves) 
potentials are compared with the experimental data from 
Ref. m. 



Ed (MeV) 


FIG. 8. (Color online) Differential cross sections ratios Rx 
for ^°Be(d,p)^^Be transfer reactions leading to ground (thick 
curves) and excited (thin curves) states of ^^Be. CX and 
CX(no subtr) predictions are given by solid and dashed- 
dotted curves, respectively. 




FIG. 9. (Color online) Differential cross section for d + ^“^Mg 
elastic (left) and inelastic (right) scattering at Ed = 70 
MeV. Results of SP, CX, and CX(no subtr) models based 
on the CH89 potential are given by dotted, solid, and 
dashed-dotted curves, respectively, while CX results based on 
the Koning-Delaroche potential are given by double-dotted- 
dashed curves. The experimental data are from Ref. [3^ . 


els show qualitatively the same behavior. The difference 
between them increases with increasing Ed and thereby 
indicates the importance of the proper fit to the two-body 
data, but does not alter the general conclusion. 

Finally I consider the d -\- ^"^Mg scattering. The 
quadrupole deformation (32 for ^^Mg determined in differ¬ 
ent experiments varies between 0.4 and 0.6 (see Ref. 
for overview). In Ref. [l^ the CX(no subtr) calculations 
were performed with (32 = 0.4 and 0.6, with the lower 
value favored by the ^■^Mg((i, d^'^Mg data and the upper 
value favored by the ^"‘Mg(fi, d'Y'^Mg* data, but no one 
of them was able to reproduce the data for both reactions 
simultaneously. Here I present in Fig. IH] the results for 
d + ^^Mg elastic and inelastic differential cross section at 
Ed = 70 MeV obtained with (32 = 0.5. The CX model 
yields the best description of the experimental data and 
fits reasonably well both elastic and inelastic cross sec¬ 
tions, while the CX(no subtr) model overestimates the 
CeX effect for the ^'^Mg((i, (i)^'^Mg but underestimates 
for the ^^Mg(d, d')^^Mg* reaction. Changing the start¬ 
ing potential from CH89 to Koning-Delaroche has only a 
small effect as demonstrated by the CX(KD) results. 


is in contrast with naive assumptions often used for the 
SF extraction by comparing DWBA-type SP calculations 
and experimental data. In Fig.|5]I confirm these findings 
of Ref. [l^ also when using an improved potential with 
CeX. For x being either CX or CX(no subtr) I show the 
ratios Rx = {da/dil)x/{da/dQ)sp in ^°Be(d,p)^^Be re¬ 
actions at 0c.m. = 0° as functions of the deuteron energy 
Ed- The deviation of Rx from the SF, which equals to 
0.855 (0.786) for the ^^Be ground (excited) state, be¬ 
comes most evident for the transfer to the ^^Be ground 
state at higher Ed- Both CX and CX(no subtr) mod¬ 


V. SUMMARY 

I considered elastic, inelastic, and transfer reactions 
in three-body nuclear systems consisting of a neutron, a 
proton, and a core. I explicitly included the core exci¬ 
tation in the exact scattering equations and solved them 
in the momentum-space framework. The calculational 
technique was taken over from Ref. [l^ but the dynamic 
input was significantly improved. 

The technical objective of this work was to develop 
a method for constructing the nucleon-core potential 
that couples ground and excited states of the core and 
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in the coupled-channel two-body Lippmann-Schwinger 
equation reproduces the predictions of the given standard 
optical potential for the elastic scattering as well as the 
DWBA predictions for the inelastic reaction. The essence 
of the method is subtracting from the original single¬ 
channel potential the explicit core excitation contribu¬ 
tions, thereby avoiding the double counting. In the usual 
case where a coupled-channel potential properly fitted to 
the experimental data is not available but the nucleus 
deformation parameters are determined in the DWBA, 
the proposed method yields the coupled-channel poten¬ 
tial including the core excitation and consistent with the 
data for both elastic and inelastic scattering. This is an 
important improvement over the simple deformation of 
the potential used previously that destroys the fit 
to the experimental data. For the p + ^°Be example it 
is demonstrated that this deviation may be significant, 
especially for the inelastic {p,p') reaction. 

The physics objective of the present work was the 
evaluation of changes in three-body observables due to 
the improvement of employed optical potentials. In 
three-body reactions involving ^°Be and ^"^Mg cores the 


most important changes were found in elastic {p,p) and 
{d,d), inelastic {d,d'), and transfer (p, d) reactions lead¬ 
ing to the core in its excited state. Transfer reac¬ 
tions (d,p) were affected less, at least at lower ener¬ 
gies. In p -I- ^^Be elastic scattering the core excitation 
effect turns out to be very small. Compared to results 
of Ref. [l^ , the improved nucleon-core potential leads to 
larger cross sections in most cases. This increases slightly 
the discrepancy between predictions and data for the 
^°Be(d,p)^^Be* transfer reaction but improves the agree¬ 
ment for the ^^Be(p, d)^°Be* reaction. In the latter case I 
also found a significant sensitivity to the parametrization 
of the standard optical potential that serves as a start¬ 
ing point in the calculations. Concerning the d + ^"^Mg 
scattering, the improved calculations are able to describe 
elastic and inelastic cross sections simultaneously with 
the same value of the deformation parameter /32 = 0.5, 
in contrast to Ref. Il9ll. Finally, I qualitatively confirmed 
the findings of Ref. [1^ that the effect of the core excita¬ 
tion in transfer reactions cannot be simply related to the 
respective spectroscopic factor. 
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